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Optical and mid- infrared neon abundance determinations 
in star-forming regions 
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ABSTRACT 

We employed observational spectroscopic data of star-forming regions compiled from 
the literature and photoionization models to analyse the neon ionic abundances ob- 
tained using both optical and mid-infrared emission- lines. Comparing Nc ++ /H + ionic 
abundances from distinct methods, we found that, in average, the abundances obtained 
via IR emission-lines are higher than those obtained via optical lines by a factor of 4. 
Photoionization models with abundance variations along the radius of the hypothet- 
ical nebula provide a possible explanation for a large part of the difference between 
ionic abundances via optical and infrared emission-lines. Ionization Correction Factor 
(ICF) for the neon is obtained from direct determinations of ionic fractions using in- 
frared emission-lines. A constant Ne/O ratio (logNe/O ps —0.70) for a large range of 
metallicity, independently of the ICF used to compute the neon total abundance is 
derived. 

Key words: galaxies: general - galaxies: evolution - galaxies: abundances - galaxies: 
formation- galaxies: ISM 



1 INTRODUCTION 

Gas phase metallicity determinations in star-forming regions 
(SFRs) have been used for the knowledge of the cosmic evo- 
lution of galaxies as well as to provide observational con- 
straints on model parameters of galaxy formation. 

In particular, the oxygen is the element most widely 
used for this purpose, because prominent emission lines 
from their main ionic stages are present in the optical 
spectra of SFRs (e.g. Hagele et al.l 120081 ; iKennicutt et al.l 
l2003l . lKrabbe fc Copettl2006f ). For the rest of elements with 
bright emission-lines (e.g. N, S, Ne, Ar) not all their ionic 
stages are observed in the optical spectrum, so their total 
abundances can only be calculated by means of Ionization 
Corre ction Factors (ICFs), proposed bv lPeimbert fc Costerd 
l|l969D . 

Among the a-elements, Ne results of great importance 



as it is one of the noble gases and it does not combine 
with itself or with other chemical species in the formation of 
molecules and, thus, dust grains. This makes the abundance 
determination of this element quite suitable for the study 
of the chemical evolution of SFRs as it does not depend on 
depletion factors onto dust grains. It is accepted that Ne 
and O are mainly produced in st ars more massive than 10 
Mq (see lWooslev fc Weaver II 19950 . Therefore, it is expected 
that a bundances of Ne and O should closely trace one to the 
other (|Crockett et al. 1120061 ). and a constant Ne/O abun- 
dance ratio over a wide range of O/H abundance must be 
found. However, recent deep spectroscopy of a large sam- 
ple of low-metallicity emission-line galaxies bv lGuseva et al.l 
(|201lT ) revealed a slight in crease of Ne/O with O/H. Similar 
results were also found by I Wang fc Liu I (|2008l ) using O and 
Ne abundances of a sample of planetary nebula e and H n 
regions. In opposite, iPerez-Montero et all (|2007l ') derived a 
new ICF for Ne using a photoionization model grid finding 
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for the high metallicity regime a slight decrease in Ne/O 
with O/H for a large sample of SFRs. 

In any case, the derivation of the total Ne abundance 
has showed to be uncertain due to the difficulties esti- 
mating its ICF. The estimation of the total Ne abun- 
dance in the optical is performed via the observation of 
the [Nelll]A3869 A emission- line, which is mainly present 
in SFRs with a high-ionization degree, and usually assum- 
ing an ICF(Ne ++ ) that is a function of the 0/0 ++ ratio 
derived also from optical emission-lines. However, the op- 
tical IKeni]_Jineis rarely observed in metal-rich objects 
(see iBresolin et al.l |2005j) , where most of the Ne is in the 
form of Ne + , so the ioni zation corrections are not well con- 
strained in this regime (|Kennicutt et al.l 120031 ). This prob- 
lem c an be alleviated by the use of infrared fine-structure 
lines (|Vermeii fc van der Hulst 1120021 ). On this context, both 
the Infrared Space Observatory (ISO, iKessler et al.l ll996T ) 
and the Spitzer Space Telescope have allowed the observa- 
tion of bright mid-infrared emission-lines in SFRs, such as 
[Neli]12.81/im and [Nem]15.56/mi, and the determination 
of the total Ne abundance simply summing the correspond- 
ing io nic fractions. For instance, IVermeii fc van der Hulst I 
1)20021) used optical and mid-IR dat a of Hn regions loca ted 
in the Magellanic Clouds (see also IVermeii et al. 1120021 ) to 
determine O, S, and Ne abundances in these objects. These 
authors found large discrepancies between Ne ++ abundances 
obtained via infrared lines with those via optical temper- 
ature determinations and pointed o ut that the classical 
appro ximation, ICF(Ne ++ )=0/0 ++ (|Peimbert fc Costerol 
Il969l), seems to underestimate the true Ne abundance (see 
also IPerez-Montero et al.ll2007l ). However, due to large scat- 
ter and uncertainty in the Ne ++ ionic fraction of their 
data, no conclusion could be drawn about the reliability 
of the optical ICF for Ne . The same results were obtained 
by iKennicutt et al.) (120031). who combined the ion ic abun- 
dances derived by Vermeil fc van der Hulst I (|2002j ) and by 
IWillner fc Nelson-Patel I (|2002T ) for H u regions located in 
M33. However, they were not able to derive a functional 
form for the ICF(Ne ++ ), due to the large scatter of the avail- 
able data and the very few measurements for low-ionization 
regions. Fortunately, a large number of mid-IR and opti- 
cal data of SFRs are currently available in the literature, 
which enables a direct determination of the optical ICF for 
Ne and a comparison among abundances obtained from dis- 
tinct methods, yielding a more reliable conclusion about the 
total Ne abundance. 

The main goal of this paper is to examine the Ne abun- 
dance of SFRs, investigating the abundance discrepancy 
found by the use of different methods and obtaining a re- 
liable ICF in the optical for this element. To do that, we 
compiled observed optical and infrared emission-line fluxes 
from the literature to derive the Ne ionic abundances. A 
photoionization model grid was also used in this analysis. 
The paper is organized as follows. In Section [5] we describe 
the observational data used along the paper and an analy- 
sis of the selected sample is performed. A description of the 
photoionization models is given in Sect. [3] The procedures 
to determine the ionic abundances are presented in Sect. [4] 
A comparison between neon ionic abundances obtained from 
distinct methods is given in Sect.[S] In Sect.[S]the ICFs com- 
puted for Ne ++ and applied to the data set are presented. 
The conclusions of the outcomes are given in Sect. [7] 
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Figure 1. Diagnostic diagrams: [Oiii]A5007/H/3 

vs. [Nn]A6584/Ha (left) and [Om]A5007/H/3 vs. 
[Sll1(A6716+A 6 731)/ Ha (right). Solid linos, taken from 
iKewlev et alj 1 120061 ), separate objects ionized by massive 
stars from the ones containing active nuclei and/or shock excited 
gas as indicated. Black squares represent HuRs and red squares 
the ELGs compiled from the literature. 



2 OBSERVATIONAL DATA 

2.1 Sample description 

Optical and infrared emission-line fluxes of a sample of 
emission-line galaxies (ELGs) - including Hn galaxies 
(HnGs), Wolf Rayet galaxies (WRGs), Blue Compact 
Galaxies (BCGs) and Dwarf Galaxies (DGs) - and Hn 
regions (HuRs) were compiled from the literature. Their 
emission-line intensities were already reddening corrected in 
the works from which we take the data. For most of the data 
this was done by using the ratio Ha/H/3 and in some cases, 
where Hcv was not present in the spectra, the reddening co- 
ef heient was obtained fr om another hydrogen lines (e.g. data 
of lCrockett et al. Il200ot ). 

The selection criterion was the presence of flux mea- 
surements of the optical [Oll]A3727, [Nem]A3869, VL/3, and 
[Om]A5007 emission-lines. We also compiled the line in- 
tensities of [Om]A4363, Ha, [Nn](A5755,A6548,A6584) and 
[Sin](A6717,A6731) when they were available. Regarding 
mid-infrared emission-lines, we selected 50 ELGs and 93 
HuRs located in the Magellanic Clouds, M101, M33 and 
M83 whose [Nen]12.81/im and [Nelll]15.56/im emission- 
lines were measured. The optical data were obtained by 
the use of long-slit, echelle and multi-object spectroscopy 
(MOS), while the infrared data were obtained with the ISO 
and Spitzer Space Telescope. 

In the optical, we used the [Om]A5007/H,S 
vs. [Nn]A6584/Ha and [O in]A5007/H/3 vs. 

[Sii](A6716+A6731)/Hq diagnostic diagrams (Figure [[]) 
to distinguish objects ionized by massive stars from 
those containing an active galactic nucleus (AGN) 
and/or gas excited by shocks. To separate the dis- 
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tinct clas ses of objec ts , we used the criteria pro- 
posed by iKewlev et al.l (|2006h . where all objects with 
log[Om]A5007/H,S < 0.61/[log[Nn]A6584/Ha -0.05] +1.3 
and log[Om]A5007/H/3 < 0.72/[log[S n]A(6717+A6731)/Ha 
-0.32] +1.3 have massive stars as their main ionization 
mechanism. From this analysis we discarded 15 objects 
yielding a sample with 522 ELGs and 212 HnRs. 

Due to the non-homogeneity of the compiled data, we 
divided the sample into five groups according to the available 
emission-lines measured for each object: 

• Group A - optical emission-lines, including the ones 
needed to directly estimate the electron temperatures and 
abundances. 

• Group B - [On]A3727, [Nem]A3869, H/3, [Om]A5007 
(the optical selection criterion) and [Om]A4363. 

• Group C - optical emission-lines except the ones sen- 
sitive to the temperature: [Oiii]A4363, [On]A7325 and 
[Nn]A5755. 

• Group D - only the [Oll]A3727, [Nein]A3869, H/3, and 
[Om]A5007 emission- lines (just fulfilling the optical selec- 
tion criterion). 

• Group E - the IR emission-lines. 

Table [l] lists the bibliographic references of the sample, the 
number of objects taken from each work, their nature, the 
observational technique used during the data acquisition, 
and the group(s) to which the data belong. The objects in 
the A, B and E groups allow the computation of the ionic 
abundances and the ICFs from several methods used along 
the paper, while the ones in the C and D groups only en- 
able a more accurate comparison between the data and the 
photoionization models to obtain a theoretical neon ICF. 

To investigate the discrepancy between the neon ionic 
abundances from optical and infrared emission-lines, we se- 
lected the objects in our sample for which we can get a 
direct estimation of, at least, an electron temperature. They 
are those belonging to the optical groups A and B, and also 
to the IR group E. The sub-sample contains 23 objects and 
they are listed in Table [2] together with the corresponding 
optical and IR references. It is worth to emphasize that only 
for the objects in this sub-sample we can estimate the oxy- 
gen and neon ionic abundances via optical lines and the neon 
ones via infrared lines. 

For the HnRs located in the Magellanic Clouds, we 
converted th e Br/3 fluxes, directly measured from the ob- 
servations bv lVermeii fc van der Hulst I &2001 ), in H/3 fluxes 
using the emissivities of IStorev fc Hummer! ( 19951 ). We did 
not use the Hucv emission line at 10.52 /mi, even though 
it is the closest hydrogen line to the mid-IR Ne ones, be- 
cause for our compiled sample it is generally blended with 
an H2 emission line. For the remaining sources we consid- 
ere d the observed H ^j fluxes. A similar procedure was done 
bv lWu et al.l (|2OO80 in their study of elemental abundances 
of BCGs using Spitzer observations. 



2.2 Analysis of the sample 

Along the years several authors have compiled spectroscopy 
data of star-forming regions from the literature in order to 
derive correl ations between macro scopic properties of these 
objects (e.g. IPilvugin et al.l 12004 ). or have suggested oxy- 
gen abundance calibrations using strong emission-lines (e.g. 




[OIII]A5007/[OII]A3727 

Figure 2. Ha/H/3 vs. [Oni]A5007/[O m]A3727. Black squares 
represent HnRs and red squares the ELGs compiled from the 
literature. All data were reddening corrected in the works from 
which we taken the data. The lines represent the theoretical val- 
ues for the low density recombination case B Balmer decrement 
Hq/H/3=3.10, 2.86 and 2. 69 for temperatu r es of 5 000, 10 000 
and 20 000 K, respectively l|Osterbrocklll989l : IStorev fc Hummer! 
ll995lV 



IPerez-Montero fc Contiiiil 120091 : iNaeao et all 120061 ). How- 
ever, the use of data sets obtained with different instrumen- 
tation and observational techniques, and including different 
objects such as whole galaxies, galaxy nuclei and individ- 
ual HnRs, i.e. non- homogeneous data, could lead to biased 
results. In what follows we investigate the possible sources 
of biases introduced in the present study due to the hetero- 
geneity of our data set. 



2.2.1 Aperture effects 

Differences between the measurement apertures used in the 
optical and IR observations could contribute significantly 
to the discrepancies f ound between optical and IR a bun- 
dance determinations (|Vermeii fc van der Hulst I I2002T ) be- 
cause many physical properties (e.g., stellar populations, 
metallicity, extinction) of galaxies vary with galactocen- 
tric radius ( Moustakas fc Kennicutti2006al ') or along nebulae 
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20001. IKewlev et all (|2005l ) presented a detailed 



analysis of the effect of aperture size on the star formation 
rate, metallicity, and reddening determinations for galaxies 
selected from the Nearby Field Galaxy Survey. They found 
that systematic and random errors from aperture effects can 
arise if fibers capture less than 20% of the galaxy light. Most 
of the SFRs in our sample can be treated as point sources 
and almost all the object extensions are observed, there- 
fore, for these objects this effect is negligible. For the spa- 
tially resolved Magellanic HnRs in our sample (see Table [2| 
this effect coul d play a more important ro l e as w as already 
pointed out by IVermeii fc van der Hulst I (J2002J). However, 
since our aim is to study the discrepancy found for the Ne ++ 
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Table 1. Bibliographic references for the compiled sample. 



ID 


Reference 


Object type 


Number 


Technique 


Group 


1 


Lee et al. ('2004') 


KISS galaxies 


13 


long-slit 


A 


2 


Vilchcz & Iglcsias-Paramo (2003) 


DG 


1 


long-slit 


A 


3 


Hagele et al. (2006) 


HnG 


3 


long-slit 


A 


4 


Hagele et al. (2008) 


HnG 


6 


long-slit 


A 


5 


Hagele et al. (2011) 


HnG 


2 


long-slit 


A 


6 


Skillman et al. (1989) 


DG 


2 


long-slit 


A 


7 


Thuan & Izotov (2005) 


BCG 


26 


long-slit 


B 


8 


Masegosa et al. (1994) 


HnG 


99 


long-slit 


A 


9 


Papaderos et al. (2008) 


HnG 


21 


long-slit 


A 


10 


Kobulnickv et al. (1999) 


HnG 


7 


long-slit 


A 


11 


Izotov et al. (2001) 


BCG 
BCG 


2 
33 


long-slit 
long-slit 


A 


12 


Izotov & Thuan (2004) a 


A 


13 


Izotov & Thuan (1998) a 


BCG 


17 


long-slit 


A 


11 


Izotov et al. (2006c) a 


HnG 


1 


long-slit 


A 


15 


Izotov et al. (1997) a 


BCG 


28 


long-slit 


A 


16 


Izotov et al. (1994) a 


BCG 
BCG 


10 
125 


long-slit 
long-slit 


A 


17 


Izotov et al. (2006b) 


A 


18 


Bergvall and Ostlin (2002) a 


BCG 


3 


long-slit 


A 


19 


Guseva et al. (2000) a 


WRG 


8 


long-slit 


A 


20 


Guseva et al. (2011) 


HnG 


59 


long-slit 


A 


21 


Guseva et al. (2007) a 


HnG 


33 


long-slit 


B 


22 


Lopez-Sanchez & Estcban (2009) 


WRG 


17 


long-slit 


A 


23 


Kobulnickv et al. (1997) 


HnR 


3 


long-slit 


A 


21 


Garnett et al. (1997) 


HnR 


7 


long-slit 


A 


25 


Kennicutt et al. (2003) 


HnR 


20 


long-slit/cchelle 


A 


26 


van Zee et al. (1998) 


HnR 


58 


long-slit 


C 


27 


Vflchez et al. (1988) 


HnR 


1 


long-slit 


A, C 


28 


Kwitter & Aller (1981) 


HnR 


6 


slit 


A, C 


20 


Bresolin et al. (2009) 


HnR 


28 


MOS 


A, C 


30 


Lee & Skillman (2004) 


HnR 


15 


long-slit 


A, C 


31 


Bresolin et al. (2004) 


HnR 
HnR 


8 

20 


MOS 
long-slit 


C 


32 


McCall et al. (1985) 


A, C, D 


33 


van Zee (2000) 


HnR 
HnR 
HnR 


2 

19 
9 


long-slit 

MOS 
long-slit 


A 


34 


Bresolin (2011) 


D 


35 


Esteban et al. (2009) 


A, C 


36 


Crockett et al. (2006) 


HnR 


10 


long-slit 


B, D 


37 


Torres-Peimbert et al. (1989) a 


HnR 


3 


long-slit 


A, C 


38 


Engelbracht et al. (2008) 


HnG 


11 


Spitzer 


E 


39 


Wu et al. (2008) 


BCG 


9 


Spitzer 


E 


10 


Gordon et al. (2008) 


HnR 


7 


Spitzer 


E 


11 


Vermeil et al. (2002) a 


HnR 


7 


ISO 


A, C, E 


12 


Lebouteiller et al. (2008) 


HnR 
HnR 


31 

21 


Spitzer 
Spitzer 


E 


13 


Rubin et al. (2008) 


E 


11 


Rubin et al. (2007) 


HnR 


24 


Spitzer 


E 



Papers from which the H/3 fluxes were obtained to compute neon ionic abundances. 



ionic abundances and this ion is not expected to be formed 
in the extended low-surface brightness zones around SFRs, 
aperture effects neither represent an important factor. 



2.2.2 Reddening correction 

There are several ways to determine the reddening correc- 
tion to be applied to the observed emission-line fluxes de- 
pending on the available information and the ad hoc assump- 
tions. For example, the use of different ratios involving hy- 
drogen emission-line fluxes, different extinction curves, and 
the use of a wide spectral range. All these facts could intro- 
duce some biases in our results. To investigate this effect, in 
Fig. [2] we plotted the Ha/H/3 ratios against the ionization 



degr ee parametrized by the |"Oih1A5007/[Ohi] A3727 ratio 
(e.g. lDors fc Copetti Il2003l ; iMorisset et al-lbool ) for the ob- 
jects belonging to our sample. In this figure we can see very 
few points (~ 3% ELGs and ~ 6% HuRs) out of the range 
of the expected values of Ha/H/3 ra t io con sidering the theo- 
retical values of lStorev fc Hummerl (|199a ) for temperatures 
between 5000 and 20000 K, typical values of star- forming re- 
gions. Thus, we can conclude that differences in the redden- 
ing correction does not affect any statistical result obtained 
in this work. 
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Table 2. Objects with neon abundances determined via optical 
and infrared emission-line intensities. 



Reference a 



Object 



Optical Infrared 



ELGs 



I Zw 18 


15 


38 


Haro 11 


18 


38 


UM420 


13 


38 


Mrk 1450 


16 


38 


NGC 4861 


15 


38 


UM 448 


13 


38 


Mrk 930 


13 


38 


Mrk 162 


13 


38 


NGC1140 


12 


38 


UM461 


13 


39 


UM462 


13 


38 


Tol2138-405 


21 


38 


IIZw 40 


19 


38 


SBS0335-052E 


11 


38 


M 101 H iiRs 




NGC5455 


37 


10 


NGC5461 


37 


10 


NGC5471 


37 


10 


Magellanic HllRs 




N160A1 


11 


11 


N160A2 


11 


11 


N157B 


11 


11 


N4A 


11 


11 


N66 


11 


11 


N81 


11 


11 



The reference number is according to the ID listed in Table [T] 



2.2.3 Instrumentation and observational techniques 

The data compiled from the literature were acquired using 
different instrumentation and observational techniques. This 
may yield an additional scatter in the abundance determina- 
tions computed along this work that might lead to biased re- 
sults. For example, if emission-lines with a broad component 
are present in low resolution spectra these could not be de- 
tected, and several physical properties obtained from these 
line s will be unreliable , such as abundance determinations 
fsee lHagele et alj|2013i . and references therein). Also, heavy 
elements abundance determinations from sp ectra with low 
signa l-to-noise ratio can be overestimated (jWesson et al. I 
I20T2T ). However, discrepancies in emission-line ratios esti- 
mated for a given object observed using different tech- 



niques ranging from about 0.03 to 0.1 dex (|Dors fc Copetti I 
I2005T ) translate into a few percent in the uncertainties of 
ionic abundance determinations. For example, if this dis- 
crepancy in the line ratio [Om](A5007+A4959)/A4363 is 
0.06 dex, the electron temperature derived using the differ- 
ent observations differs by about 400 K, yielding a log(Ne/0) 
variation of about 0.1 dex. This value is of the order of 
the intrinsic scatter found for uniform samples studied us- 
ing data obtained wit h the same instrumentation (see e.g. 
iKennicutt et al.l I2003T ). When no systematic effects in the 
calibrations and data reduction exist, the physical proper- 
ties derived using data obtained with different observational 



techniques should be the same. For example, the physical 
parameters estimated for the star forming knots of Haro 15 
are in very go od agreement when derived using long-slit and 
echel le data (jHagele et al.l |2012| ; lLopez-Sanchez fc Estebanl 
120091) . 



2.2.4 Nature of the sources 

Our sample is composed by objects with different morphol- 
ogy and sizes, nevertheless they all have young and massive 
stellar clusters as their main ionization sources. The shape 
of the emission-line spectra of the objects in our sample is 
dominated by the ionizing fluxes from these massive stars. 
Therefore the same measurement and analysis techniques 
can be applied to derive the temperatures, densities and 
chem ical composition of their interstellar gas (jHagele et al.l 
l2006h . 

In general, the HllRs and ELGs are located in the 
same zone in diagnostic diagrams (see e.g. I Baldwin et al.l 
Il98lf ) implying that similar ionization mechanisms are tak- 
ing place in these two different kind of objects (see Fig. [lj|. 
Other ionizing mechanisms could be present according to 
the evolutionary stages and the nature of t he objects, suc h 
as galactic winds due to supernovae (e.g. iTang al.l 120071 ). 
contribution to the hydrogen emission fluxes by diffuse ion- 
ized gas (e.g.lOev et al.ll2007l ; iMoustakas fc Kennicuttll2006l ; 



I Calzetti et ifeooA . a ncTToT the presence of obscured I o¥s 
( Kewlev et al.l 120011 ). or even the effects of the presence of 
multiple kine matical components in the emission-line pro- 
files (see e.g lHagele et all |2010| ; iFirpo et al.l |2010| . l201ll ; 
lAmorin et alj 120121 . and references therein). The statistical 
contribution of these mechanisms is similar for the differ- 
ent kind of objects and they do not yield any bias on the 
physical conditions of the objects obtained in the present 
work. 



3 PHOTOIONIZATION MODELS 

To estimate the ICF for Ne via photoionization models, 
we compared the results of a grid of models performed 
using the Cloudy/8.0 code (JFerland et al.lll998h with the 
observational data sample. We used data of the groups A, 
B, C, D, and E which include all objects with optical and IR 
emission-lines available (i.e. 734 objects). The diagnostic di- 
agrams [Oiii]A5007/[Oii]A3727 vs. [Nn]A6584/[0 n]A3727, 
[Nem]A3869/H/3 vs. [On]A3727/H/3, 

[Nii]A6584/Hq vs. [Sn]A6720/Ha, and 

[Nem]15.56/i/[Neii]12.81/im vs. [Oin]A5007/[Oii]A3727 
were used for the comparison (see Fig. [3]). The ICF val- 
ues were obtained from the model results which better 
reproduce the observational data. 

The model grid was bu ilt following the same pro- 
cedures as iDors et al. I (|2011f ). with metallicities of Z = 
2.0,1.0,0.6,0.4,0.02,0.01 Z e , and logarithm of the ioniza- 
tion parameter U from —3.5 to —1.5 dex with a step of 
—0.5 dex. For each model, the ionizing source was as- 
sumed to be a stellar cluster with a spectral energy dis- 
tribution obtained from the stellar population synthesis 
code Starburst99 (jLeitherer al.lll999f ). with an upper mass 
limit of 100 Mq and an age of 1 Myr. For the hypotheti- 
cal nebula we adopted a constant electron density N e — 200 
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cm 3 , plane-parallel geometry, and solar abundance ratio 



log(Ne/O) = -0.61 (oxygen from JAllende Prieto et alJliooH 
and Ne from iGrevesse fc Sauva]Hl998T ). In the models with 
Z = 2 Zq, a value of log(Ne/0)=— 1.0 was assumed, in or- 
der to fit the data in this hi gh metallicity regime. This v alue 
is about the same found by IPerez-Montero et al.l (|2007l ) for 
very high metallicities. A complet e description of the pho- 
toionization models is presented in lDors et al. I (|201lh . 

In Fig. |3] the diagnostic diagrams containing the ob- 
served emission-line intensities and the photoionization 
model results are shown. The typical errors of the obser- 
vational emission-line ratios are about 10 per cent. It can be 
seen that the majority of the observational data falls within 
the regions occupied by the models when all the considered 
emission lines are in the optical range (see both upper and 
lower left panels). However, in lower right panel of Fig. [JH 
can be seen that the models produce lower values of the in- 
tensity line-ratio [Ne h i] 15.56/jm/ [Ne h] 12.81/im. This p rob- 
le m was pointed out bv | Perez-Montero fc Vflchea (J2009J) and 
by iMorisset et al.1 (|2004T ) . They found similar results from 
photoionization model grids built assuming different atmo- 
sphere models as photoionization sources. They concluded 
that this effect could not be due to uncertainties in stel- 
lar atmosphere models and must be due to problems in the 
atomic data of the mid-IR lines. 



4 DETERMINATION OF IONIC 
ABUNDANCES 

We considered two methods to determine the Ne ++ , Ne + , 
+ , and ++ ionic abundances in our sample. The method 
that uses direct estimations of the electron temperatures 
(this method will be called the Visible-lines method) and 
the one that uses infrared emission-lines (this method will 
be called the IR-lines method). 



4.1 Visible-lines method 

For those objects with the appropriate optical emission-line 
measurements, we calculated the electron temperature (T e ) 
from the observed line-intensity ratio 7?o3=[Oin](A4959 + 
A5007)/A4363 for the high ionization zone (refereed as 
£3). To do that w e used the fitting function derived by 
lHagele et al.l (|2008l ) based on the temden routine of the 
nebular package of the iRAljj software: 



£3 = 0.8254 - 0.00024157?o3 + 



47.77 
R03 ' 



(1) 



with £ in units of 10 4 K. The electron density (N e ) was com- 
puted from the ratio of [S n]A6716/A6731 using the temden 
routine and £3 values given by the expression above. 

The ++ and Ne ++ abundances were computed using 
£3 and following the relations given bv lPerez-Montero et al.l 
(|2007l ): 



12 + log( 



(V 



H + 



log( 



1(4959) + 1(5007) - 
J(H/?) 



+ 6.144 



1 Image Reduction and Analysis Facility, distributed by NOAO, 
operated by AURA, Inc., under agreement with NSF. 
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+— 0.55 log £3 



(2) 



and 



(3) 



+ 1^_ 0.504 log i 3 . 
£3 

The + ionic abundances were computed assuming 
N c = 100 cm" 3 and an electron temperature of the low ion- 
ization zone (refereed as £2) derived from the theoretical 
relation: 



tz 1 = 0.693 fcT 1 +0.281 



(4) 



obtained from the models described by 

IPerez-Montero fc Diazl (|2003l ). Measurements of the 
[Oll]A7325 auroral emission line, sensitive to the [On] 
electron temperature, are available for only 10-15 per- 
cent of our sample. Since we used the observational 
sample in a statistical way, we preferred to derive 
all the temperatures following the same procedure. 
Besides, almost all objects present high or moderate 
excitation degrees [(/[O ni]A4959+A5007)/(/[O ll]A3727+ 
7[Q H I1A4959+A5007) > 0.4, as was defined by iPilvuginl 
|2000j. Therefore the impact of + /H + variations to the 
total oxygen abundance and to the ++ /(0 ++ +0 + ) ratio 
(used in the neon ICF determinations) is small (~18%) in 
any case, and from a statistical point of view negligible. 
Only in those cases where the [Om]A4363 and [On]A7325 
auroral emission-lines were not available, we computed the 
+ ionic abundances assuming the electron temperature 
directly estimated from the observed line-intensity ratio 
i?N2=[Nn](A6548 + A6584)/A5755: 



£2 = 0.537 - 0.000253/i N 2 + 



42.13 

i?N2 



(5) 



fitting function given bv lHagele et al.1 (|2008f h For few objects 
(< 1%) £2 was derived using this equation. 

To derive the + ionic abundances we used the expres- 
sion given by Hagele and collaborators: 

,0+, 



12 + log( 



H+ 



log (&)+ ^92 



+ - 



/(H/3) 

1.583 



£2 



0.681 log £ 2 + log(l + 2.3n c ),(6) 



where n e — N c /(W 4 cm~ 3 ). As pointed out by these authors, 
the expressions given above are valid for temperatures from 
7 000 to 23 000 K. The + ionic abundances calculated using 
£2 values from Eqs. [4] and [5] differ by about 0.2 dex. The 
data set used to determine the ionic abundances using this 
method belongs to groups A and B and is comprised by 
579 objects (~ 80% of the sample). All + , ++ and Ne ++ 
ionic abundances obtained via Visible method used along 
this study were computed from this data set. 



4.2 IR-lines method 

The Ne + and Ne ++ ionic fractions can be determined us- 
ing the intensities of the [Nell]12.81^m and [Nelll]15.56^m 
emission-lines, respectively, and following a simil ar method- 
ology to the one presented bv lForster et al.l (|200lh . The ionic 
abundance of an element can be given by 
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Figure 3. Diagnostic diagrams containing observational data taken from the literature (all groups) and the results of the grid of 
photoionization models described in this work. The solid lines connect curves of iso-(0/H), while the dashed lines connect curves of 
iso-U. The colour code for the values of the oxygen abundanc es (tracer of the metallicity ) are the same in all panels and are indicated 
in the upper left panel. The solar oxygen abundance refers to lAllende Prieto et al.l 1 120011 ) and corresponds to 12+log(0/H)=8.69. The 
values of log U are indicated in each panel. Points are the observational data. The typical errors (bars not shown) of the observational 
emission-line ratio are about 10 per cent. 



n x+i _ h(X +i )N e j x{H+) 

71 H + /a(H+) Ne j A (X + i)' 



(7) 



where n x +i and n H + are the densities of the ions X +1 
and H + , 7a(X +1 ) is intensity of a given emission- line emit- 
ted by X +1 , I\m+) i s the intensity of a reference hy- 
drogen line, while jAfx+M an d Ja(h+) are the emissivity 
values given by the IONIC routine of the nebular pack- 
age of the ir af software w h ich u se s the Ne atomic pa - 
rameters from iBadnefl et all (J2006t). iGriffin, et all l|200ll) . 



Kaufman fe Sugar I (| 19861 ), iButler fc Zeippen I (| 19941) , and 
Mendoza I |l983|). These emissivity values are assumed to 
be constant in all abundance determinations because these 



vary less than 5 % over a large temperature range (see the 
pioneer work bv lSimpsonlll975l ). 

Using this method any error in the determination of 
these emissivities directly translates in a systematic shift in 
the derived Ne + and Ne ++ ionic abundances. Taking into 
account all these assumptions we obtained 



Nc+ 

h+ 

and 



/(12.81/xm) 



x 1.322 10" 



Ne ++ _ /(15.56pm) 



H+ 



7(H/J) 



x 6.323 10" 



(8) 



(9) 
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12 + log(Ne ++ /H + ; 



Visible— lines method 



Figure 4. Comparison between ionic abundance of the Ne + ~'~/H~'~ 
derived using Visible-lines and IR-lines methods. Points represent 
estimations for the objects presented in Table [2] where squares 
are estimations using H/3 fluxes and triangles the ones using 
Br/3 fluxes. Solid line represents equality of the two estimates. 
Dashed line is the equality line shifted by the average of the dif- 
ferences between the ionic abundances. Typical errors in optical 
and infrared estimations are of a bout 15% and 10%, respectively 
l|Vermeii fc van der Hulst 1120021 ). 




[Nelll]15.55/irri 
[Nelll]A3869S 



^dust/^van Hoof 

Figure 5. Logarithm of the luminosity predicted by our models 
for some hydrogen and neon emission-lines versus the dust in- 
ternal abundance in relation to the default value assumed in the 
Cl oudy code, i.e. l og O/ H fa —4. The grain model is described 
in Ivan Hoof et al.l I I2001I) . Points are results for models consider- 
ing different dust abundances connect by the curves. The results 
for different emission-lines are indicated. The Z dust /Z v!ln Hoof = 
represents the models free of internal dust. 



5 COMPARISON OF NEON IONIC 
DETERMINATIONS 

With the aim to compare the Ne ++ /H + ionic abundances 
derived using the Visible-lines method and those using the 
IR-lines method, we plotted in Fig. [4] the results for the 23 
objects listed in Table [2] the only ones for which we are 
able to apply both methods. These objects are those that 
simultaneously belong to groups A or B, and E. We can note 
a large discrepancy and scatter. Comparing the Ne ++ /H + 
ionic abundance determinations obtained via the Visible- 
lines method with those based on the IR-lines method, we 
found discrepancies of about a factor of 4 in average, be- 
ing Ne ++ /H + underestimated when using the Visible- lines 
method (see Fig. [4]). 

The origin of this discrepanc y is uncertain and several 
explan ations have been proposed. I Vermeil fc van der Hulst I 
(2002) compared optical and infrared ionic fractions ob- 
tained for a sample of HnRs located in the Magellanic 
Clouds. These authors found considerable discrepancies be- 
tween determinations via these two methods, which were 
attributed to the difference in the area of the sky covered by 
the infrared and optical observations. However, as discussed 
above, aperture effects are not probably the origin of this 
huge discrepancy regarding a high-excitation ion as Ne ++ . 

Extinction might noticeably affect the comparison 
between the abundances obtained from IR and optical 
emission-lines. The blue optical [Nem] emission-line is more 
absorbed by dust than the mid-IR ones. Hence, if the nebu- 
lar emission is highly reddened by dust the optical line does 
not trace all the Ne in the inner parts of the SFRs, yield- 



0.4 

C(H/3) 



Figure 6. Differences between Ne +_l " ionic abundances obtained 
via the IR-lines and the Visible-lines methods as a function of the 
nebular reddening constant C(H/3) for the majority of the objects 
listed in Table [2] 
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ing an optical ionic abundance lower than that derived from 
the mid-IR lines. In the same way as for the Ne emission- 
lines, this effect could be enhanced by the use of H/3 or Br/3 
emission lines instead of Hua to derive ionic abundances in 
the mid-IR. In order to verify the extinction influence on 
the emission-lines, models with different internal dust abun- 
dance were built. In these, the only effect of grains on the 
continuous opacity was considered, i.e. the gas heating by 
grain emission was not considered in the calculations. The 
dust abundance was linearly scaled with the default value 



stellar atmospheres. Therefore, this could be symptomatic 
of inaccurate atomic data involving these emission-lines and 
can be on the basis of the found discrepancies when de- 
riving abundances from optical and mid-IR emission lines. 
In particular, the atomic parameter which can affect the 
abundance determinations is the collision strength for the 
[Nem]A15.56^m emission-line. However, it is not probable 
that uncertainties in this parameter be the main cause of 
the discrepancy found in our study. Along the years differ- 
ences in neon collision strength have been found. For ex- 



assumed in the Cloud y code, i.e. log O/H ~ —4. The grain a mple, the collis ion strength for [Ne ill] Pi — >■ P% derived 



model is described in Ivan Hoof et al.l (J200l|). In Fig. [5] the 
predicted luminosity values (L/Lfrecdust) for some hydrogen 
and neon emission-lines in relation to the ones predicted 
by the model free of dust (Z& ust /Z va , n Hoof=0) are shown. 
We can see that H/3, Br/3 and Hua have about the same 
decreasing with the dust abundance, reaching ~4% maxi- 
mum difference in their luminosities for Zd us t/-Z vall Hoof=10. 
A similar behaviour is found for the neon emission-lines 
with a difference of about 6%. Even though these differ- 
ences are within the observational measurement errors, we 
have used these predicted emission-line intensities to calcu- 
late the Ne ++ abundance following the same procedure pre- 
sented in Section [3] The largest difference found in Ne ++ 
abundances derived from IR and the Visible methods was 
0.3 dex. Therefore, this result shows that, at least, the ex- 
tinction caused by the internal dust is not the main source 
of the neon ionic discrepancy. 

To verify if there is an observational correlation between 
the extinction and the abundance discrepancy for the 23 
objects listed in Table [2l we plotted in Fig. [6] the differ- 
ence (D) between the Ne ++ ionic abundances obtained via 
the IR-lines and the Visible-lines methods versus the nebu- 
lar reddening constant [C(H/3)] compiled from the literature 
and used by the authors to correct their observational data. 
There is no evident correlation between the discrepancy and 
the extinction found from the Balmer decrement. 

For five objects of our sample: N81, N4A, N66, N160A1, 
and N160A2, both Br/3 and H/3 fluxe s were obtained directl y 
from spectroscopy observa tions by Vermeil et al. (120021) ■ 



iDufour fc Harlow I (|l977al lbl) and iHevdari-Malaveri et al.l 



dl998l.l2002r). We use t hese data and the [Ne m]15.55/im from 
IVermeii et al. I (|2002h to calculate Ne ++ /H + ionic abun- 
dances of these objects using both hydrogen lines. We found 
neon ionic estimations in agreement by about 0.2 dex, in- 
dicating that dust absorption is not the main source of the 
discrepancy, at least for these objects. Another observational 
support t o this result arises from the data presented by 
IWu et al.l (|2008l ). These authors pointed out that most of 
the galaxies observed by them, which include 50% of the 
ELGs presented in Table [2] do not show a strong 9.7/^m 
silicate feature and, consequently, have low dust extinction. 
Moreover, as noted in Fig. [4l for the objects whose both 
Br/3 and H/3 fluxes were measured, the ionic difference is 
also found, indicating that the differential absorption of op- 
tical and infrared emission lines is not the main cause of the 
discrepancy. Thus, we can assume that this effect is negligi- 
ble for these objects. 

As can be seen in Sect. [3l and was already not i ced b y 
iPerez-Montero fc Vflcliezl l|2009l ) and iMorisset et all (|2004t) , 
the photoionization models can not reproduce diagnostic di- 
agrams based on Ne mid-IR lines independently of the used 



by iBlaha I (|l969l) is f2=0.581. The fl valu es derived in the 
IRON project bv lButler fc Zeippen I |l994h range from 0.481 
to 0.778 for temp eratures ranging from 10 3 to 10 5 K. Later 
computations by iMcLaughlin fc Belli (|2000r i yield a lower 
SI valu e by abou t a 20% than the one pr oposed by IBlaha I 
(|l969j). Recently. IMcLaughlin et al.l (|2011f ) have recomputed 
the collision strengths for some neon lines using a small 
scale 56-level Breit-Pauli calculation and a large-scale 554- 
levels R-matrix Intermediate Coupling Frame Transforma- 
tion (ICFT). They showed that different Q values are ob- 
tained by the use of these suppositions, in the sense that for 
56-level Breit-Pauli calculation Q ranges from about 0.5 to 
1.0 and for the 554-levels ICFT ranges from about 0.4 to 0.8. 
It must be noted that these results are in good agreement 
with earlier derivations. Nevertheless, e ven assuming that 
the collision strength for the [Nelll]A15.56/im emission-line 
can vary by a factor of two, it is not enough to conciliate 
the ionic abundance discrepancy. 

Finally, another supposition is the presence of electron 
temperature fluctuation amplitude inside the ionized gas. 
It can lead to an underestimation of the ionic abundances 
obtained via optical collisionally-excited lines when these 
fluctuations are not considered since these lines ar e strongly 
dependent on the electron temperat ure (see e.g. IPeimberj 
ll967i . l2003l ; lKrabbe fc Copettill2005h . This problem can be 
alleviated using abundances based on IR emission-lines, 
whose fine-structure transitions h ave a weak dependence 
on the electron temperature (e.g. iDinerstein et al. I Il985l; 



Rubin et af] l200ll ; iNollemberg et all I2002J ; iGarnett et all 



2004T ). Since the electron temperatures are related with the 



chemical composition of the emitting gas, temperature fluc- 
tuations are essen tially equivalent to abund ance variations 
across the nebula (jKingdon fc Ferlandll 19981 1 . 

To test the effect of the presence of temperature fluctu- 
ations and/or abundance variations, we built a grid of pho- 
toionization models using Cloudy. These models are similar 
to the ones described in Sect. [jJ] and presented in in Fig. [3] 
but considering the presence of chemical inhomogeneities in 
the theoretical nebula and spherical geometry. The plane- 
parallel geometry assumed in the later models can not be 
used to test variations along the radius of the hypothetical 
nebulae, since makes the inner radius much larger than the 
thickness of the cloud and the results about independent on 
electron temperature fluctuation amplitude. In this case we 
considered a nebula with Z = Zq, with an inner radius of 
4pc and an outer radius where the temperature falls below 
4000 K (~30pc), and a number of hydrogen-ionizing pho- 
tons emitted by the ionizing star cluster of 10 50 s . It was 
assumed that the metallicity varies along the radius of the 
nebula as a sine wave, with period P ranging from 0.01 to 
2pc, including 0.01, 0.1, 0.25, 0.5 and 1 pc. Also a model 



10 Dors et al. 




IR— lines method 

Visible — lines method 



0.5 



1.5 



Period (pc) 



Figure 7. Results of photoionization models considering abun- 
dance variations along the radius of the nebula. Ne++/H+ ionic 
abundances obtained via the Visible-lines and the IR-lines meth- 
ods (dashed- and solid-lines, respectively) as a function of the 
period P of the abundance variations. 
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Figure 8. Temperature structure of the theoretical nebula from 
our photoionization models. For each model, R/R a is the radius 
normalized to the outermost one (Ho). Bottom panel: The curves 
represent result models with different values of P (period of the 
sine wave which the mctallicity varies along the radius of the 
nebula) as indicated. Top Panel: The curves represent model re- 
sults with no abundance variations but with different distances 
between the central ionizing source and the inner face of the gas, 
repres ented by d, as indicated. The t 2 as defined by IPeimberti 
119671 ) and the electron temperature £3 for each models arc shown. 



with no variations and solar abundance was considered. We 
considered a set of models with an abundance amplitude of 
0.8 dex, which takes metallici t y valu es between 2 x Zq and 
0.3 x Z Q . iKingdon fc Ferlandl (|1995l ) built grids of photoion- 
ization models using similar sine wave variations but applied 
to the total hydrogen density across the nebula. 

We used the emission-line intensities given by our mod- 
els for each value of P to calculate the Ne ++ /H + ionic abun- 
dances applying the Visible-lines and the IR-lines methods 
following the same procedure as described in Section [4] 



Thus, we obtain two estimations of Ne ++ /H + for each 
model. Fig. [7] shows these estimated values of Ne ++ /H + 
plotted as a function of the period P of the abundance vari- 
ations. As it can be seen when no variations are considered, 
the Ne ++ /H + ionic abundances estimated by each method 
differ by less than 0.1 dex, which is within the errors of the 
methods. When abundance variations are present the differ- 
ences could reach about 0.5 dex, being the ionic abundances 
via the IR-lines method the largest one. This value is very 
similar to the one found in Fig. [4] between the estimated 
Ne ionic abundances derived applying the Visible-lines and 
the IR-lines methods to the observational data. Moreover, 
these differences, as expected, depend on the metallicity am- 
plitude of the models, increasing (decreasing) with an in- 
crement (decrement) in the amplitude (models not shown). 
Therefore, the temperature fluctuations studied through the 
abundance variations across the nebula could be a probable 
explanation for the discrepancy between the estimated Ne 
abundances by these two methods. In Fig.[8](bottom panel), 
we have plotted the temperature structure [j of our models 
with P=2, 0.01 and no fluctuations versus the radius nor- 
malized to the outermost one. In this figure we also show, 
for each P, the values of the temperatur e fluctuation pa- 
rameter (i 2 ) as defined by IPeimberti (|1967l ) and the electron 
temperature £3 using the predicted values of R03 and the 
Eq.[T] 

We found, as expected, that abundance variations 
translates into temperature fluctuations. As can be seen in 
Fig. |HJ £3 decreases when larger values of P are assumed, 
however, the average temperature weighted by the electron 
density over t he volume of th e nebula <T> (not shown), 
as defined by IPeimberti (|1967l ). is about constant with P. 
The model with P=0.01 pc shows a constant temperature 
fluctuation amplitude of about 9000 K, while the ampli- 
tude showed by the model with P=2 decreases with the ra- 
dius. Interestingly, the t 2 value, when no fluctuations are 
con sidered, is about the sa me found in som e HllRs (see 
e.g. iKrabbe fc Copettll2002l ) and HiiGs (e.g. lHagele et all 
2006). On the other hand, the t 2 values that explain the dif- 
ferences found when comparing ionic abundances obtained 
via recombination and collisional excited emis sion-lines (see 
e.g. IPeimberti 120031 : IGarcia-Roias et al.l 120061) are in good 
agreement with those given by the models with P differ- 
ent than zero. A more detailed study about spatial electron 
temperature fluctuations i n ionized nebu lae from photoion- 
ization models is given bv lCopettil (|2006l) . 

Along decades observations have failed in found a 
direct evid ence of chemi c al in h omogeneities in H 11 re- 
gions (e.g. lOliveira et al.1 [20081: iLebouteiller et~al| 120081 : 
IGarcia-Roias et all ^OrJoT Krabbe fc Copettil 12003 ) and 
in planetary nebulae (e.g. Krabbe fc Copettil 120051 ; 
iRubin et aT]|2002l : iLiu et alj|2000» . Thus, an another source 
of temperature fluct uations can be the presence of density 
variations in the gas. iRubinl |l989j) showed that large-scale 
variations in electron density in the gas do arise electron 



2 In Cloudy code the nebula is divided in concentric zones with 
physical conditions across them about constant. This causes the 
zone number in each computed model to vary. The electron tem- 
peratures of these zones are shown in Fig. \E\ See Cloudy manual 
for more details. 



Ne abundance determinations in star-forming regions 11 



temperat ure flu c tuatio ns. Considering the Hn region 
models of iRubinl (|1985T ). this author pointed out that if N c 
varies by a factor of about 10, the electron temperatur e 
can range up to 2000 K (see als o lKingdon fc Ferlandlll995h . 
Recently, IMesa-Delgado et al.l ( 2012T ) presented results 
from integral field spectroscopy of a region near to the 
Trapezium Cluster and found the presence of high-density 



(N c 



4 x 10 cm ), which are much higher than 



the typical values ob tained either for ga lactic and extra- 
galactic regions (see ICopetti et al.l 120001 ). Unfortunately, 
a larger number o f obse rvations, such as the ones by 
IMesa-Delgado et al.l (|2012i ). are needed to confirm the high 
scale density fluctuations in ionized nebulae. 

Another likely reason for the presence temperature fluc- 
tuations are different proximities of gas to the ionizing O 
stars, since these are mainly responsible for the ionization 
structure. To test this, we built models considering different 
distances, d, between the central ionizing source and the in- 
ner face of the gas, where d ranged from 0.1 to 7 pc. These 
values are consistent with the ones used bv lSimpson et al.l 
(|1997T ) in a study of the Galactic Hn region GO. 18-0.04. 
Again, the predicted emission-lines were used to calculate 
the Ne ++ /H + ionic abundances applying the Visible-lines 
and the IR-lines methods, following the same procedure as 
described in Section [4] We found a 0.0-0.2 dex difference in 
Ne ++ /H + . In Fig. [8] (top panel) the temperature structure of 
the models is shown. We can note a hot region very near the 
central star, however, this region is too small to affect the 
derived abundances, since no significant changes among the 
models neither in the t 2 nor in £3 values are obtained. Mod- 
els with a more complicated ge ometry, with many s tars dis- 
tributed through the gas (e.g. Ercolano et al.ll2007l ). might 
cause the temperature variations and consequently help to 
explain Ne ++ /H + ionic discrepancies. 



6 NEON ICFS 

6.1 ICFs determinations 

The ICFs stand for the unseen ionization stages of each el- 
ement. For the Ne ++ it is defined as 



ICF(Ne ++ ) = Ne/Ne"* 



(10) 



The ICF for Ne ++ can be directly computed using only 
[Nell]12.81/im and [Nelll]15.56/im emission-line fluxes. As- 
suming that the total Ne abundance is 



Ne _ Ne + Ne ++ 

h" ~ "fF + h+ ' 

and using Eq. [10] we obtained 
ICF(Ne ++ ) = l + ^L. 
Using Eqs. [5] and [5] we found 



(11) 



(12) 



F 



[NcII]A12.81/im 



Fn 



(13) 



ICF(Ne ++ ) = 1 + 2.10 

" , [NeIII]A15.56Mm 

Therefore to estimate the Ne ICF using the IR lines it is not 
necessary the use of hydrogen emission-line fluxes. Thus, we 
were able to use all the objects belonging to group E (143 
objects) . A histogram containing the ICF values lower than 
10 is shown in Fig. [9] It can be seen that almost all the ICF 
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Figure 9. Histogram containing the values of the ICF(Ne++) 
obtained using the IR emission-line intensities and Eq. ll3l for the 
objects belonging the group E and with ICF values lower than 
10. 



values (~ 85%) are lower than 5. Very high values (> 10) are 
also found. The objects that present the largest ICF values 
are IC342 («257), NGC5236 and NGC 5253 (~58), Searle5 
(»55), NGC 3628 (-45), and NGC 2903 («32). Also high 
ICFs were obtained from dat a of some H nRs in M 83 ob- 
served by iRubin et al.l (J2007T ) . These values (not shown) 
could reflect the presence of cool ionizing sources in these 
objects. 

A classical way to determine the neon ICF takes advan- 
tage of its relation with the ionic oxygen ++ /(0 + + ++ ) 
ratio. From a fit in our photoionized model results presented 
in Sect. [3] we obtained: 



ICF(Ne ++ ) modc i 



0.741 - 0.08x 



0.393 



(14) 



where x=0 ++ /(0 + + ++ ). This relation is shown in 
Fig. [10] In this figure are also shown the ICFs obtai n ed us - 
ing pho toionization models b y IPerez-Montero et al.l (|2007l ) 
and by llzotov et al.l (|2006bT ) . the classical approximation 
[ICF(Ne ++ )=0/0 ++ ] which the ionic abundances were de- 
rived from Visible-lines method and using our data sample, 
and the ICFs obtained from the IR-lines method. 

Fitting the results obtained using the ICF via IR- 
lines method and x via Visible-lines method, represented 
by points in Fig. 1101 we obtained the following expression: 



ICF(Ne ++ )ia = 2.382 - 1.301x + — , 

x 



(15) 



valid for x Js 0.4. Unfortunately, only for 23 objects 
in our sample we have measurements of the Ne infrared 
emission-line intensities and we were able to determine the 
++ /(0 + +0 ++ ) ratio via the Visible-lines method (see Fig. 
I10[) . It must be noted that this sample comprehends objects 
in a large range of metallicity (7.1 < 12 + log(0/H) < 8.50 
dex) as well as a large range of ionization degree (0.4 < 
++ /(0+ + ++ ) < 1.0). 
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Photoionization models 
Perez-Montero et al. (2007) 
Izotov et al. (2006) high Z 
zotov et al. (2006) inter, and 

Ne ++ /Ne = ++ /0 Classicc 
R — lines method 




x = ++ /(0 ++ + + ) 

Figure 10. Relation between ICF(Ne++) and x=0++/(0+ + 
++ ) ionic abundance ratio. Solid lines correspond to our fit- 
tings to the ICF for Ne and those found by different authors as 
indicated. Points represent the data whose ICFs were determined 
via IR-lines method and listed in Table [2] 



6.2 Application of the Neon ICFs 

To test the differences in Ne/H total abundance yielded by 
the use of different ICFs, we used the observational data 
to compute Ne ++ /H + ionic abundances via the Visible- 
lines method and applied different ICFs to obtain the 
Ne total abundances. Fig. [TT] shows the comparisons be- 
tween the Ne/H derived using the ICF from the IR- 
lines method (Eq. I15[) an d those values estimat e d app lying 
the ICFs de r ived by : (a) IPerez-Montero et al.l J2007J 1. (b) 
llzotov et al.l (|2006bf ). (c) photoionization models (Eq. I14p 
described in this work, and (d) the classical approxima- 
tion [ICF(Ne++) = 0/0++]. It can be seen that, in 
general, the expression from the IR-lines method gives 
higher ICF values, and consequently higher Ne/H values, 
by about 0.08-0.13 dex than those obtained using ICFs from 
IPerez-Montero et all J2007f ). lizotov et all (|2006bl ). and from 
our photoionization models . For low-e xcitation objects, we 



find that the high metallicity models bv llzotov et al.l (|2006bf ) 
are in agreement with our theoretical relation. Finally, the 
classical ICF approximation yields Ne/H values in conso- 
nance with the ICFs via the IR-lines method. 



6.3 Ne/O vs. O/H 

A very important issue 
between the Ne/O ratio 
dance , for which there is 
l|2008l ). using oxygen and 
of planetary nebulae and 
ratio increases with O/H 
similar resu l t was found 
llzotov et all |2006rl ) for a 



is the study of the relation 
and the total oxygen abun- 

not a consensus. IWang fe Liu I 

neon abundances of a sample 
HuRs, found that the Ne/O 
in both types of nebul ae. A 
by iGuseva et all (|201lh and 

large sample of low- metallicity 



ELGs. IPerez-Montero et al.l (|2007l ). using a photoioniza- 
tion model grid to derive the Ne ICF found, for a large 
sample of SFRs, that the Ne/O behaviour agrees with 
the assumption of a constant value for the low metallicity 
regime (12+log(0/H)<8.2) but shows a slightly decrease for 
the high metallicity regime (12 +Iog(Q/H)>8.2). Moreover, 
IWillner fc Nelson-PateJl (|2002l '). using mid-IR derivations of 
the neon abundances of a sample of HnRs located in M33, 
obtained a flat neon abundance gradient as a function of the 
galactocentric radius, as well as a decrement of the Ne/O ra- 
tio when the O/H increases. 

To investigate this issue, we applied the different ICFs 
for Ne shown in Fig. [TO] to the ionic Ne ++ abundances es- 
timated using the Visible-lines method to obtain total Ne 
abundances for the compiled sample. We then combined 
these values with the total oxygen abundances to estimate 
Ne/O ratio. Figs. 1 121 and 1 131 show these values as a function 
of 12+log(0/H) for our ICFs and those from the literature, 
respectively. In these figures we also show the logarithm of 
the average Ne/O value derived using the different ICFs, and 
their standard deviations. These average values are consis- 
tent in all the cases with the adopted solar one. From a 
visual inspection of the distributions of our observational 
sample in these diagrams, can be seen that there is a fiat 
relation between the logarithms of the Ne/O and the O/H 
ratios for the whole studied metallicity range, from about 7 
to 8.5 dex. We also perform a fitting to these data assuming 
a linear regression with slope a and the regression constant 
6, and without taking into account the individual errors. In 
Table[3]we list the coefficients of the fittings. We obtain null, 
positive and negative slopes for the different ICFs. However, 
in all the cases these slopes are close to zero. In the worst 
case, a — —0.08, a metallicity variation equal to the complete 
metallicity range (1.5 dex, equivalent to a factor of about 
32) implies a variation of the Ne/O ratio equal to 0.12 dex 
(equivalent to a factor of about 1.3). This Ne/O variation 
is very similar to the standard deviation of the data. We 
also analysed the statistical dependence of Ne/O with O/H 
using the Spearman's rank correlation on the data, find- 
ing that the Ne/O is statistically constant independently 
of the metallicity value. Another test was done consider- 
ing different metallicity regimes. We computed the average 
Ne/O value for the high (12+log(0/H)>8.2) and the low 
(12+log(0/H)<8.2) metallicity regime for each considered 
ICF. In Table [3] these results are listed. None variation of 
the Ne/O with O/H is noted for these different regimes and 
for each ICF. Therefore, we can conclude that the Ne/O and 
the metallicity have a flat relation for this sample. This flat 
Ne/O vs. O/H behaviour for the whole range of O/H anal- 
ysed is a robust test of nucleosynthesis prediction and sup- 
ports a very limited (if any) oxygen depletion in dust. The 
oxygen abundances i n grains are poorly known. For example, 
llzotov et alj (|2006bf ) interpreted that the slight increase of 
Ne/O with the metallic ity is due to depletion of abo ut 20% 
of O onto grains. Also, IPeimbert fc Peimbertl (|201Cl ) found 
an increasing depletion of O atoms of about 0.1 dex with 
increasing O/H for three objects in a large range of metal- 
licity. From our data we derived a dispersion of Ne/O for a 
fixed O/H value of ~0.10 dex. Therefore, none conclusion 
can be obtained about the oxygen depletion on dust. 

On the other hand, it must be noted that the scat- 
ter of the Ne/O increases with the metallicity, showing the 
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Figure 11. Comparison between total Ne/H abundances obtained from the Ne++/H+ ionic ones using the Visible-lines method and 
applying different ICFs as indicated. Points represent estimations for the objects belonging to the groups A and B. The top panel of 
each diagram shows the difference between the estimations using the methods considered. The average value D of this difference and the 
dispersion are shown in each plot. Solid lines represent equality of the two estimates. 



distribution of the da ta a triangular-l i ke sh ape. This re- 
sult was also found bv lKennicutt et al.l (120031 ) for HnRs in 
M101, NGC2403 and dwarf irregular galaxies. They argued 
that this effect can reflect high sensitivity to the local ra- 
diation field as Ne ++ and ++ become minor constituents 
of the nebular material when the + /0 increases, which 
is equivalent to a metallicity increment. As the metallicity 
of the gas material increases the electron temperature de- 
creases, thus, the Ne ++ and ++ emission-line intensities 
becomes weak. Therefore, the increment of scatter going to 
the higher abundances should be a normal behaviour admit- 



ting that errors in the measurements of the sensitive lines 
are found to be larger (statistically) for higher abundance 
HnRs. In fact the triangular shape mentioned is probably 
just a product of this. We have also tested the dependence 
of Ne/O on the temperature fluctuations, since neon and 
oxygen have slightly different temperature coefficients. Us- 
ing the earlier result models and computing Ne/H and O/H 
following the same procedure described in Sect. [4] we found 
that log(Ne/0) is -0.53 and -0.36 for P=0.0 and 0.01, re- 
spectively. Thus, temperature fluctuations would not be a 
source of systematic errors in the Ne/O derivation. 
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Table 3. Coefficients of the linear regressions log(Ne/0)=a X [12+log(0/H)]+b fitted to data plotted in Figs,IT2land ll3l Columns 2 and 
3 list the slope and the constant coefficient, respectively, together with their formal errors. Columns 4 and 5 are the mean values of the 
log(Ne/0) for the high and low metallicity regimes considered. 



log(<Ne/Q>) 



ICF(Ne++) 



high Z 



low Z 



IR-lines method 
Photoionization models 
Izotov et al. (2006b') 
Perez-Montero 

ICF(Ne++)=0/0++ 



+0.02 ±0.01 
-0.02 ±0.01 
±0.06 ±0.01 
-0.08 ±0.01 
±0.04 ±0.01 



-0.77 ±0.11 
-0.53 ±0.10 
-1.17 ±0.12 
-0.08 ±0.14 
-1.06 ±0.12 



-0.63 ±0.11 
-0.74 ±0.13 
-0.70 ±0.12 
-0.83 ±0.13 
-0.64 ±0.11 



-0.62 ±0.08 
-0.70 ±0.09 
-0.75 ±0.10 
-0.75 ±0.13 
-0.65 ±0.09 



o 



R — lines method 
log(<Ne/0>) = -0.62 ±0.09 



Photoionization models 
log(<Ne/0>) = -0.71 +0.10 
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Figure 12. Relation between log(Ne/0) and O/H using the 
Visible-lines method and our two different ICF(Ne+ + ) as in- 
dicated. Points represent estimations for the objects belong- 
ing to the groups A and B. The solid line represents the 
Ne/O abundance rat i o der ived using the oxyg en abundance from 
Allcn de Prieto et al.l IJ20011 ) and the neon from lGrevesse &: Sauvall 
1119981), while the das hed line represents the value derived by 
iDrake fc Testa I 1120051 ) using Chandra X-ray spectra of a sample 
of nearby solar-like stars. The average and the standard devia- 
tion of log(Ne/0) taking into account low and high metallicity 
regimes are shown in each plot. 



Finally, we tested the Ne/O dependence on the ioniza- 
tion degree x=0 ++ /(0 + + ++ ). To do that, we used the 
compiled sample to obtain the Ne/O abundances via Visible- 
lines method and the x values. The ICF for Ne from IR-lines 
method (Eq. I15|) was considered. In Fig. [14] (bottom panel) 
we can see that higher Ne/O values are found for objects 
with high excitation. This result is independent of the ICF 
considered. Although the neon ICF and O/H ratio decrease 
with x, the behaviour of the Ne/O ratio is due to the incre- 
ment in the Ne ++ /H + with x, as can be seen in Fig. 1141 (top 
panel). 



Izotov et al. (2006J 
og(<Ne/0>) = -0.73 ±0.12 



■ ■ i 



Perez-Montero et al. (2007) 
log(<Ne/0>) = -0.77 ±0.13 



■- ^^r-r^M^ ^i ^:- 



ICF(Ne + + ) = 0/0 + + 
og(<Ne/0>) = -0.65 ±0.10 



7.5 



8.5 



12±log(0/H) 



Figure 13. As Fig.[l2]for the ICF(Ne++) from the literature as 
indicated. 



7 CONCLUSIONS 

In this work we show that the ionic fractions Ne ++ /H + ob- 
tained via optical direct detections of the electron temper- 
ature are underestimated by about a factor of 4 in relation 
to those estimated via mid-infrared emission-lines. We dis- 
cussed several possible causes responsible of this discrep- 
ancy, including: 

(i) Different apertures between optical and 
mid-IR observation^ as was already suggested by 



IVermeii fc van der Hulstl (|2002t ). This effect should 
not be of great importance for the most part of the anal- 
ysed sample, which behaves as point-like sources. Besides, 
this effect has minimum effect in a high-excitation ion such 
as Ne ++ . 

(ii) Differential extinction effects. Mid-IR emission-lines 
are possibly tracing a deeper ionization structure than opti- 
cal ones. Besides, due to no reliable measurements of the 
Hucv hydrogen recombination line at 10.52 fim since it is 
blended with an Eb line, ionic abundances were derived in 
relation to H/3 or Br/3 emission-lines. We did not find any 
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Figure 14. Bottom panel: Relation between log(Ne/0) and 
the ionization degree x=0 ++ /(0 + + ++ ) obtained using the 
Visible-lines method, the compiled sample and the ICF for Ne 
from the IR-lines method. Points represent estimations for the 
objects belonging to the groups A and B. Top panel: Relation 
between estimations of 12+log(Ne ++ )/H + ) and x. 



correlation between the degree of discrepancy and the inner 
extinction found from the Balmer decrement. 

(iii) Uncertainties in the atomic data related with the 
infrared emission of both [Ne n] and [Ne ill] . These emission- 
lines present some problems as it is evidenced by inspect- 
ing the right lower panel of Fig. 1 in which observational 
data and photoionization models are shown. According to it 
the models predict systematically lower values of the mid- 
IR [Ne ill] /[Nell] emission-line ra tio than observations. This 
discrepancy was already found b v lMorisset et al.l (J200J) and 
IPerez-Montero fe Vilchej (|2009l ) and apparently it is not re- 
lated with the shape of the used spectral energy distribution 
in the models. The values for the [Nem]A15.56/im emission- 
line collision strength found in the literature along decades 
differ by about 50%. Even when we assume that collision 
strength for the [Nem]A15.56/im varies by a factor of two, 
it is not enough to conciliate the neon ionic abundance dis- 
crepancy. 

(iv) Variations of chemical abundances across the neb- 
ula. Using a grid of photoionization models we showed that 
fluctuations with an amplitude of 0.8 dex in the total abun- 
dance (for both O and Ne) can produce the degree of dis- 
crepancy between the ionic Ne ++ abundances derived from 
optical and mid-IR emission-lines. 

Using ionic Ne abundances estimations from mid-IR 
emission-lines, we obtained an expression for the ICF of the 
Ne ++ as a function of the ++ /(0++0 ++ ) ratio. These 
ICFs obtained through the use of empirical and theoretical 
derivations of Ne ++ from mid-IR lines come from a quotient 
of emissivities and therefore minimize the impact of points 



(ii) and (iii) mentioned above, which are more critical for 
the abundance in absolute value. We then employed a large 
sample of observational spectroscopic data of star-forming 
regions compiled from the literature to analyse the depen- 
dence of the Ne/O abundance ratio with O/H. We found 
that Ne/O is about constant with O/H (log Ne/O w -0.70) 
for the whole metallicity range considered. This result is in- 
dependent of the approximation adopted for the ICF, and 
in all the cases the average Ne/O estimated for our sample 
is consistent with the adopted solar value. 
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